Complete fusion of 9 Be with spherical targets 
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The complete fusion of 9 Be with 144 Sm and 208 Pb targets is calculated in the coupled-channels 
approach. The calculations include couplings between the 3/2~, 5/2~, and 7/2~ states in the 
K — 3/2 ground state rotational band of 9 Be. It is shown that the B(E2) values for the excitation 
of these states are accurately described in the rotor model. The interaction of the strongly deformed 
9 Be nucleus with a spherical target is calculated using the double-folding technique and the effective 
M3Y interaction, which is supplemented with a repulsive term that is adjusted to optimize the 
fit to the data for the 144 Sm target. The complete fusion is described by in-going-wave boundary 
conditions. The decay of the unbound excited states in 9 Be is considered explicitly in the calculations 
by using complex excitation energies. The model gives an excellent account of the complete fusion 
(CF) data for 9 Be+ 144 Sm, and the cross sections for the decay of the excited states are in surprisingly 
good agreement with the incomplete fusion (ICF) data. Similar calculations for 9 Be+ 208 Pb explain 
the total fusion data at high energies but fail to explain the CF data, which are suppressed by 20%, 
and the calculated cross section for the decay of excited states is a factor of three smaller than the 
ICF data at high energies. Possible reasons for these discrepancies are discussed. 

PACS numbers: 24.10. Eq,25.60.Pj, 25. 70.-Z 



I. INTRODUCTION 

The influence of breakup on the complete and incom- 
plete fusion of weakly bound nuclei with stable targets 
is currently being studied at many experimental facilities 
around the world. Experiments with unstable nuclei are 
particularly challenging because of weak beam currents 
and poor statistics. Fortunately, there are several light 
elements that are both stable and weakly bound and they 
provide the opportunity to study the influence of breakup 
on fusion with good statistics. 

A good example of a stable and weakly bound nu- 
cleus is 9 Be and its fusion has been measured with 208 Pb 
and 144 Sm targets [H-i]. The simplest view of 9 Be is a 
strongly deformed three-body system consisting of two a 
particles held together by a weakly bound neutron. The 
Q-value for the a + a + n three-body breakup is -1.574 
MeV. In both experiments it was possible to separate the 
complete from the incomplete fusion. 

One of the most dominant features in coupled-channels 
calculations of the fusion of the strongly deformed 9 Be 
nucleus with a stable target is the excitation of states in 
the K=3/2 ground state rotational band of 9 Be. The 
large quadrupole deformation of 9 Be (derived from the 
measured quadrupole moment of the ground state) im- 
plies that conventional calculations that are based on a 
deformed Woods-Saxon potential may become unrealis- 
tic, for example, if the curvature corrections to the ion- 
ion potential 0] , which are due to the deformation of the 
reacting nuclei, are ignored. These problems are over- 
come in the following by using the double-folding tech- 
nique [1, @ to calculate the Coulomb plus nuclear inter- 
action between the deformed 9 Be nucleus and a spherical 
target. 

Another interesting feature of 9 Be in the non-zero spin 
of the ground state. This feature was pointed out in Ref. 



[H, and it was recommended that the spin of the 3/2~ 
ground state and the 5/2~ and 7/2~ excited states of 
9 Be should be treated explicitly in the calculations. In 
particular, the fusion cross section should be calculated 
for each of the initial magnetic quantum numbers, m = 
±1/2 and ±3/2, of the 3/2~ ground state of 9 Be, and 
the average cross section should be compared to measure- 
ments. In the rotor model one can easily calculate the 
necessary matrix elements from the multipole expansion 
of the total interaction between projectile and target. 

The structure and the parametrization of the one-body 
density of 9 Be is discussed in the next section. The calcu- 
lation of the double-folding interaction between 9 Be and 
a spherical target is presented in Sect. III. A model for 
calculating the complete fusion cross section is presented 
in Sect. IV. The results of coupled-channels calculations 
of the fusion of 9 Be± 144 Sm and 9 Be± 208 Pb are presented 
in Sect. V, and Sect. VI contains the conclusions. 



II. STRUCTURE OF BE 

The nucleus 9 Be behaves like an almost perfect rotor 
with respect to quadrupole excitations of the K = 3/2 
ground state rotational band with spins I" = 3/2 _ , 5/2~, 
and 7/2~ . This can be seen by comparing the measured, 
reduced transition probabilities for quadrupole transi- 
tions to the results obtained from the expression, Eq. 
(4.68a) of Ref. 0], 



B(E2,KI ->. KI') = (IK 2Q\I'K) 2 , 



(1) 



which applies to a perfect rotor. Here Qq is the intrinsic 
quadrupole moment. The measured quadrupole moment 
of the 3/2" ground state of 9 Be is 5.29(4) fm 2 which 
translates into the intrinsic quadrupole moment Qq = 



2 



TABLE I: The measured reduced transition probabilities 
B(E2, 3/2" ->■ I~) for exciting the K=3/2 ground state ro- 
tational band of 9 Be [ij are compared to predictions of the 
rotor model. Also shown are the excitation energies Ej and 
decay widths Ti [9j. 
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3/2 _ 








14.0(2) 


14 


5/2 _ 


2.429 


0.78 


40.7(6) 


36 


7/2" 


6.38 


1210 


18.9(3) 


20 


sum 






73.6(11) 


70 



TABLE II: Density parameters for 9 Be, 144 Sm, and 208 Pb. 
The measured RMS charge radii [llj shown in the last col- 
umn are reproduced, and so is the intrinsic quadrupole mo- 
ment of 9 Be, Qo = 26.45(20) fm 2 0. The last line shows 
the adopted neutron (y) density parameters for 208 Pb; the 
'RMS-exp' radius is estimated from Skyrme HF calculations 




Nucleus 


R (fm) 


a (fm) 


P 


< r 2 > x / 2 


RMS-exp 


9 Be 


2.08 


0.375 


1.183 


2.540 


2.52(1) 


144 Sm 


5.829 


0.54 





4.941 


4.947(9) 


208 Pb(7T) 


6.60 


0.546 





5.500 


5.503(2) 


20S Pb(u) 


6.82 


0.546 





5.66 


5.66 



26.45(20) fm 2 , according to Eq. (4-69) of Ref. 0. In- 
serting this value into Eq. (JXJ) one obtains the reduced 
transition probabilities that are shown in the last column 
of Table I. They are in good agreement with the measured 
values shown in the 4th column of Table I. The sum of 
the B- values is shown in the last row of Table I. The sum 
is 5QQe 2 /(16"7r) in the rotor model and differs from the 
experimental value by about 5%. 

A great advantage of the rotor model is that it can 
be applied to calculate the transition matrix elements 
in cases where they have not been measured, for exam- 
ple, for the 5/2" to 7/2~ transition, and the quadrupole 
moments of the excited states. One can also calculate 
matrix elements of the total interaction between 9 Be and 
a spherical target from a multipole expansion of this in- 
teraction as discussed in Sect. III. 

The excitation energies and widths of the three states 
are shown in the 2nd and 3rd columns. It is noted that 
the width of the 7/2" state is very large which implies 
that the state, if excited, may decay by particle emis- 
sion before the fusion with the target takes place. This 
possibility is investigated in Sect. IV and V. 

There are other low-lying states in 9 Be but they are 
not expected to play any significant role in heavy-ion 
collisions. Thus the spin-orbit partners of the ground 
state rotational band, i. e., the 1/2", 3/2~, and 5/2~ 
states, will be ignored because spin-excitations are weak. 
The excitation of positive parity states, starting with the 
lowest l/2 + state, is also weak and will be ignored. 



A. Density parametrization 

The densities of the deformed projectile and spherical 
target nuclei will be parametrized by the expression 

, _ r l + cosh(fl/q) 
P[T ' ] ~ ° cosh(r/a) + cosh(i?/a) ' ( ' 

where C is a normalization constant, R is the radius, 
and a is the diffuseness. The radius of the deformed, ax- 
ially symmetric projectile depends on the direction with 



respect to the symmetry axis. It is parametrized as 

R{9') = Ro[l + /3 2 Y 20 (9')], (3) 

where 9' is the angle between the position vector r and 
the direction e of the symmetry axis. 

The advantages of the parametrization @ are that 
it is similar to a Fermi function at larger values of r 
and it is well behaved as a function of 9' for r — > 0, 
where it approaches an orientation independent constant. 
Another advantage of Eq. (ffl ) is the analytic properties 
it has for spherical nuclei [1(|. For example, the Fourier 
transform is an analytic function, and the expression for 
the root-mean-square (RMS) radius, 

^) = \{R 2 + \{a*) 2 ), (4) 

is an exact relation (see the appendix of Ref. [13].) These 
features will be utilized for the spherical target nuclei, 
144 Sm and 208 Pb. The parameters that have been cho- 
sen are shown in Table II. They were adjusted so that 
the measured RMS charge radii were reproduced. 
The parameters for the neutron (i>) densities were as- 
sumed to be the same as for protons (tt), except in the 
case of 208 Pb, where a slightly larger radius is used to 
accommodate for the neutron skin of this nucleus. The 
adopted skin thickness: S np =< r 2 >\l 2 — < r 2 >^ 2 
~ 0.16(6) fm was chosen because it falls in the midst of 
values predicted by Skyrme Hartree-Fock (HF) calcula- 
tions [l2| • Moreover, it is also consistent with the skin 
thickness S np = 0.16 ±(0.02) stat ± (0.04) sys fm that has 
been extracted from antiprotonic 208 Pb atoms [11]. The 
parameters for 9 Be are determined below. 



B. Multipole expansion of density 

The density of the deformed nucleus is expanded on 
Legendre polynomials, 

pM') = EM r m( cos ( 0, ))< (5) 
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FIG. 1: (Color online) Correlation between the intrinsic 
quadrupole moment and RMS charge radius of 9 Be. The 
curves were obtained by varying @2 for the fixed radius param- 
eter R — 2.08 and the three values of the diffuseness indicated 
in the figure. 



where p\{r) is calculated numerically, 



P\{r) = 



2A + 1 



dx P\(x)p(r). 



(6) 



The multipole expansion of the Fourier transform of the 
density is 



p(k) = JVVa(*) P A(C0S(^), 



(7) 



where 8' k is the angle between k and the direction e of 
the symmetry axis, and 

/>oo 

px(k) =4ir dr r 2 p x {r) j\{kr). (8) 
Jo 

The above expressions are used in the next section to 
calculate the double-folding potential. They are also used 
to calculate the electric multipole moments of the de- 
formed charge density p c (r, #'), 



M(E\p) = M(EX) Y Xfl (e), 



(9) 



where 



M(EX) 



4tt 



2A + 1 



dr r A+2 p cX {r). (10) 



The intrinsic quadrupole moment Qq is traditionally de- 
fined as 2M(E2). 



C. Calibration of the density of 9 Be 

The parameters of the density of the deformed 9 Be nu- 
cleus were adjusted so that both the intrinsic quadrupole 
moment and the RMS charge radius agree with experi- 
ments. That was achieved as follows. The mean square 
charge radius of 9 Be, 



(r) = T 



dr r 4 p c ,A=o(r), 



(11) 
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FIG. 2: (Color online) The shape of 9 Be, Eq. Q, derived by 
calibrating the density to reproduce the measured intrinsic 
quadrupole moment Qo and RMS charge radius. 



and the intrinsic quadrupole moment, 



Qo = 2M(E2) 



47T 

T 



dr r 4 PcA=2(r). (12) 



were calculated as functions of the deformation parame- 
ter /?2 for a fixed radius R — 2.08 fm and for three values 
of the diffuseness. The results are shown in Fig. Q] as a 
correlation between the RMS charge radius and Qq. It 
is seen that the curve which is based on the diffuseness 
a = 0.375 fm passes trough the experimental values, and 
agreement with both values is achieved for /3 2 = 1.183. 
This is the value that will be used in the following, and 
the shape it produces according to Eq. Olooks almost like 
two touching a particles as illustrated in Fig. [2j In fact, 
the intrinsic quadrupole moment of 9 Be, Qo = 26.45(20) 
fm , is almost identical to the calculated quadrupole mo- 
ment of the unbound nucleus 8 Be. The published value 
obtained in Variational Monte Carlo (VMC) calculations 
is 26.6(3) fm 2 [lij]. The intrinsic quadrupole moment of 
10 Be, on the other hand, is slightly smaller; the value one 
obtains from the measured B{E2) value of the lowest 2 + 
excitation @ is Qq = 22.9 fm 2 . 



III. DOUBLE-FOLDING POTENTIAL 

Having adopted the rotor model for 9 Be and deter- 
mined the densities of the projectile and the spherical 
144 Sm and 208 Pb targets, one can now use the double- 
folding technique to calculate the potential that will be 
used in the coupled-channels calculations. The double- 
folding potential is defined by 



U(v) 



J dn J dr 2 p{ri,e) p T (r 2 ) u(|r 2 +r-n|), 



(13) 



where v is the effective nucleon-nucleon interaction and 
r is the relative distance between projectile and tar- 
get. The target density px is assumed to be spherical 
whereas the density of the projectile 9 Be is deformed and 
parametrized as described in the previous section. 



4 



The double-folding potential is calculated most conve- 
niently from the Fourier transforms of the densities ac- 
cording to the expression 



U(r) = 



dk 



p(k) pr(-k) v(k) e 4kr . (14) 



Inserting the expression ([7]) for the deformed projectile 
and spherical target densities into Eq. (fT4)) one obtains 

U(r) = U(r,e') = J> A (r) *M«»(0')), (15) 



A. Matrix elements 

Having expressed the total interaction U(r, 6') 
(Coulomb + Nuclear) in terms of the multipole expan- 
sion Eq. (p~5|) . one can now easily calculate the diagonal 
as well as off-diagonal couplings between states in the 
ground state rotational band of 9 Be. All one needs to 
calculate is matrix elements of the Legendre polynomi- 
als, 

P A (co S (0')) = ££$(r) ^o(e). 



where 6' is the angle between r and e and 

U x (r) = ^ J dk k 2 p x {k) PT (k) v(k) j x (kr). (16) 

The double-folding calculation of the ion-ion poten- 
tial and its multipole expansion, Eq. (|15[) . will be based 
on the M3Y effective interaction, supplemented with a 
repulsive term that simulates the effect of nuclear in- 
compressibility. This method has been applied previously 
by Mi§icu and Greiner [6j to calculate the fusion between 
spherical and deformed nuclei. It was also used in Ref. 
[1 51 ] to explain the hindrance in the fusion of spherical 
nuclei at extreme subbarrier energies. 

The repulsive term in the effective TV TV interaction is 
parametrized as a contact interaction, 



V NN'( r ) = VrS(r), 



(17) 



and the densities that are used in the associated double- 
folding calculation have the same radius as shown in Ta- 
ble I but the diffuseness a r is usually chosen much smaller 
[lf|. The value chosen here is a r — 0.3 fm. 

The Coulomb interaction interaction can also be gen- 
erated from Eq. ([TO simply by replacing v(k) by the 
Fourier transform Aire 2 /k 2 of the proton-proton Coulomb 
interaction, and replacing the nuclear densities with the 
charge densities p c of projectile and target. The results 
has the same form as Eq. (jl"5j) . 



L/ c M')=5>A C (r) Pa(cos(0')), 



(18) 



where 



1 f AlTP 2 

U?(r) = ^ dk k 2 p c ^(k) Pc , T {k) — j x (kr). (19) 



The matrix elements between different states are 



(KI'M'\P x {cos(6'))\KIM) 



21 
21 



7^ Y. D lUr)(IMXp\I'M') (IKX0\I'K). (21) 



The calculation is even simpler in the rotating frame ap- 
proximation, which is used in the coupled-channels cal- 
culations described in the next section. In this approxi- 
mation one assumes that r (the relative distance between 
projectile and target) defines the z-axis. The angle 9' is 
then identical to the angle e of the symmetry axis with 
respect to the z-axis. Since D^ (z) = S^q this implies 
that p=0 is the only non-zero term in Eq. (|2"Tj) . 

The total potentials one obtains for the two systems 
9 Be+ 144 Sm and 9 Be+ 208 Pb are shown in Fig. O Shown 
are for each system the monopole potential (solid line) 
and the entrance channel potentials for the magnetic 
quantum numbers m=l/2 and 3/2, of the 3/2~, K=3/2 
ground state of 9 Be. All three potentials were obtained 
with the strength v r = 410 MeV fm 3 of the repulsive 
effective NN interaction (which is determined in subsec- 
tion V.A.) 

The magnetic quantum numbers m=l/2 and 3/2 used 
in Fig. [3] refer to a z-axis that points in the direction of 
the relative position of projectile and target. The to=3/2 
channel therefore corresponds to an orientation where 
the tip of the deformed 9 Be points towards the target, 
whereas the m=l/2 corresponds to the belly pointing to- 
wards the target. Consequently, the Coulomb barrier for 
the to=3/2 entrance channel is lower than the barrier 
for m—1/2. Another observation is that the potential 
pocket is deeper for m=l/2 than for m=3/2. This is a 
consequence of a larger radius of curvature and a stronger 
nuclear attraction for the m—1/2 belly configuration. 



For large separations of projectile and target this inter- 
action approaches the usual monopole-multipole interac- 
tion, 



c Z T e 2 M(E\) 



(20) 



IV. MODEL OF 9 BE INDUCED FUSION 

The cross sections for the complete fusion of 9 Be with 
a heavy target are calculated in the coupled-channels ap- 
proach. The complete fusion is simulated by ingoing wave 
boundary conditions that are imposed in all channels at 
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FIG. 3: (Color online) Entrance channel potentials for 
9 Be+ 144 Sm and 9 Be+ 208 Pb, respectively, obtained for the re- 
pulsive strength v r — 410 MeV fm 3 . The solid curve is the 
monopole potential; the dashed curves are the entrance chan- 
nel potentials for the magnetic quantum numbers m — 1/2 
and 3/2. The ground state energies of the two compound 
nuclei, 153 Dy and 217 Rn, are indicated. 



the minimum of the pocket in the entrance channel po- 
tential. The coupled equations are solved in the rotat- 
ing frame approximation [l7T - ll"9| , where the z-axis points 
in the direction of the separation vector r between the 
reacting nuclei. This approximation is also called the 
isocentrifugal approximation fi~7j because the centrifugal 
potential is the same in all channels and equal to the 
centrifugal potential in the entrance channel. The total 
magnetic quantum number to is also preserved in this ap- 
proximation (see Ref. [13] )• Since the ground state of the 
target is a + state and the ground state of 9 Be is a 3/2~ 
state, one would have to solve the coupled equations four 
times, for m = ±1/2 and ±3/2. For symmetry reasons it 
is actually sufficient to calculate the fusion cross section 
twice, for to = 1/2 and 3/2, and the average fusion cross 
section will be compared with data. 

The effect of the decay of the excited state in 9 Be is 
included explicitly in the coupled equations by employing 
the complex excitation energies E] — ijTj. The decay 
widths Tj are shown in the 3rd column of Table I; it 
is seen that the width of the 7/2~ state is very large. 



I 




FIG. 4: (Color online) Measured cross sections for the com- 
plete fusion of 9 Be+ 144 Sm || are compared to calculations for 
m = 1/2 and 3/2, and the average cross section (solid curve). 
The strength of the repulsive interaction was set to v r = 410 
MeV fm 3 . 



The dominant decay mode of the 7/2~ state is neutron 
emission, and 55% of it populates the 2 + excited state of 
8 Be Q. The latter state has an excitation energy of 3.03 
MeV and a large width of 1.51 MeV, with an exclusive 
decay into two a particles ||. The two a particles are 
emitted back to back in the 8 Be rest frame, so if one of 
them is emitted towards the target, the other a partner 
will recoil away from the target nucleus and will most 
likely escape. It is therefore assumed that the decay of 
the 7/2~ state will not lead to complete fusion (CF), 
and the CF will be calculated as described above from 
the ingoing flux, whereas the decay is registered in the 
absorption cross section. 

The data for the complete fusion of 9 Be± 144 Sm Q are 
compared in Fig. [3] to the results of coupled-channels 
calculations. The calculated CF cross sections were de- 
rived from the ingoing flux as described above. The top 
dashed curve is the cross section for the to = 3/2 mag- 
netic substate, the lower dotted curve is for to = 1/2. It 
is seen that the curve for to = 3/2 dominates the CF at 
all energies, consistent with the lower Coulomb barrier 
for this magnetic substate (see Fig. The data should 
be compared to the solid curve which is the average of the 
CF cross sections for m = 1/2 and 3/2. The comparison 
is discussed in more detail in subsection V.A. 



A. Incomplete fusion 

The decay of the 7/2~ state will end up in the breakup 
of 9 Be. The precise outcome of the decay in terms of 
incomplete fusion or breakup is not so clear. It would 
require a multi-cluster description to follow the two a 
particles after the decay. As mentioned earlier, it is un- 
likely that both a particles fuse with the target because 
they are emitted back to back. However, it is possible 
that one of them will fuse with the target nucleus and 



6 



1 



S 



1000 
100 
10 

1 

0.1 

1000 
800 
600 
400 
200 








(A) > 




A * 




/ ¥ 


% 144 C 

Be on Sm 


i .4 






CF — 
ICF • 



25 30 35 40 45 



TF > 






. CF ► 


— ♦— 




ICF > 
















J* (B) 









25 



35 
(MeV) 



40 



FIG. 5: (Color online) Measured complete (CF) and incom- 
plete (ICF) fusion cross sections for 9 Be+ 144 Sm 0] are com- 
pared in (A) to the calculated cross sections for CF (solid) 
and absorption (dashed curve). The linear plot in (B) also 
shows the total fusion (TF) cross section. 



lead to incomplete fusion (ICF). It should also be em- 
phasized that there are other sources of ICF, for example, 
the neutron transfer from 9 Be, which are not included in 
the coupled-channels calculations presented here. 

The calculated cross section for the decay of the excited 
states will in the following be referred to as the absorption 
cross section. In view of the above discussion one should 
not expect that the absorption would account for the 
measured ICF cross sections but it is clearly of interest 
to compare the two cross sections. The experimental 
total fusion (TF) cross section is the sum of the CF and 
ICF cross sections. It will be compared to the calculated 
TF cross section, which is the sum of the CF and the 
absorption cross sections. 



V. COMPARISON TO MEASUREMENTS 

The results of coupled-channels calculations that are 
based on the model presented in the previous sections 
are compared in the following to the data for the com- 
plete (CF) and incomplete fusion (ICF) of 9 Be+ 208 Pb 
[H, [| and 9 Be+ 144 Sm Q. Both targets are spherical, 
closed shell nuclei and the excitation of these nuclei is 
relatively weak compared to the excitation of 9 Be and 
they will be ignored. The fusion with the 144 Sm target 
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FIG. 6: (Color online) The \ 2 P er data point for the com- 
plete (CF) and total fusion (TF) of 9 Be+ 144 Sm are shown as 
functions of the strength v r of the repulsive interaction, Eq. 
(Ell- 



is discussed first because the couplings are weaker in this 
case and the adopted model is therefore expected to be 
more successful. This case will also provides the oppor- 
tunity to calibrate the repulsive part of the effective NN 
interaction. 



A. Fusion of 9 Be+ 144 Sm. 

The calculated cross sections for the fusion of 9 Be with 
144 Sm are compared in Fig. [5j\. to the data of Ref. [![. 
The measured and calculated cross sections for complete 
fusion (CF) are seen to be in good agreement. This was 
achieved by adjusting the repulsive part of the effective 
NN interaction which is used in the calculation of the 
double-folding nuclear potential. The best fit to the data 
is obtained for the strength v r ss 410 MeV fm 3 , and that 
is the value that will be used in the following. The x 2 
per data point is shown in Fig. [S] as a function of the 
strength of the repulsive interaction. There is another 
solution with a small x 2 f° r v r < 300 MeV fm 3 but it is 
unphysical because it produces a pocket for 9 Be+ 208 Pb 
that is deeper than the energy of the compound nucleus. 

The dashed curve in Fig. (5j\ shows the calculated 
cross sections for the decay of the excited states of 9 Bc; 
it is in surprisingly good agreement with the ICF data. 
The good agreement may be accidental but it could also 
indicate that the decay of the excited states of 9 Be is 
the main source of ICF for the 144 Sm target. The best 
fit to the data for total fusion (TF) is also achieved for 
v r w 410 MeV fm 3 . This consistency of the CF and TF 
is illustrated in Fig. [5] in terms of the x 2 P er data point. 

The measured and calculated fusion cross sections are 
compared in a linear plot in Fig. [5)3. It is seen that 
the different components of the measured and calculated 
fusion cross sections are in good agreement. In particu- 
lar, the CF is suppressed by about 10% compared to the 
total fusion (TF) at high energies, both experimentally 
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FIG. 7: (Color online) Measured cross sections for the com- 
plete (CF) and incomplete fusion (ICF) of 9 Be+ 208 Pb [| are 
compared to the calculated cross sections for CF (solid curve) 
and absorption (dashed curve). The linear plot in (B) also 
shows the total fusion (TF) cross sections. 



and in the calculations. This suppression is caused in the 
coupled-channels calculations by the decay of the excited 
states of 9 Be. Without any decay in the coupled-channels 
calculations, the fusion cross section obtained from the 
ingoing-wave boundary conditions would be close to the 
measured TF cross section. In other words, the suppres- 
sion of the CF compared to the TF requires some sort 
of absorption mechanism, and the decay mechanism sug- 
gested here seems to provide a natural explanation. Let 
us now investigate whether this mechanism can explain 
the data for the Pb target. 



B. Fusion of 9 Be+ 208 



Pb. 



The results of the coupled-channels calculations of the 
fusion of 9 Bc+ 208 Pb arc compared in Fig. [7]to the data of 
Ref. |2[. The calculations are similar to those presented 
above for 9 Be+ 144 Sm. It is seen in Fig. 03 that the 
calculated absorption cross section (due to the decay of 
excited states) can only explain 1 /3 of the measured ICF 
cross section at high energies. The suppression of the 
CF compared to the TF cross section is about 30% in 
the experiment @, whereas the calculations only show 
a 10% suppression. There are evidently other sources of 
ICF in collisions of 9 Be with a 208 Pb target, besides the 



decay of excited states considered here. 

A 30% suppression of the CF data was observed in 
Ref. by comparing to coupled-channels calculations. 
The calculations were based on a deformed Woods-Saxon 
potential but did not consider the effects of incomplete 
fusion nor the decay of excited states. It was shown that 
a scaling of the calculated fusion cross section by a factor 
of 0.7 leads to a very good agreement with the CF data 
at all energies. In Fig. [7] it is sufficient to multiply the 
CF calculation by a factor 0.8 in order to match the CF 
data at high energies. The reason for the smaller scaling 
factor is that the decay of the excited states has already 
taken care of a 10% reduction. 

It is often necessary to employ a weak, short-ranged 
imaginary potential in order to be able to reproduce the 
fusion data of stable nuclei at high energies by coupled- 
channels calculations. This is particularly the case when 
the calculations are based on a shallow entrance chan- 
nel potential [lOj. Since the potentials shown in Fig. 
[3] are relatively shallow it is of interest to see what is 
the effect of a weak imaginary potential on the fusion of 
9 Bc+ 208 Pb. Let us therefore choose the potential 



W{r) = Wo [1 + exp((r - R w )/a w )} 



(22) 



with a w = 0.2 fm and R w = 9.5 fm so that it acts near the 
minimum of the potential pockets shown in Fig. [3] The 
strength Wo was adjusted to optimize the fit to the CF 
data. The best fit is shown in Fig. [8j^; it was achieved 
for W = -2.5 MeV and has a X 2 /N=1A. 

The absorption cross section shown in Fig. [5J due to 
the combined effect of the imaginary potential, Eq. (|22l) , 
and the decay of the excited states, is in good agreement 
with the ICF data at high energies but the discrepancy 
is large at low energies. The discrepancy indicates that 
the breakup leading to ICF must take place at larger 
separations of projectile and target than assumed in the 
potential, Eq. (|2"2")l . It may be possible to construct a 
more realistic imaginary potential (of the volume plus 
surface type) but that idea will not be pursued here. 

There are other reaction channels that could be a 
source of ICF. Examples are the dissociation of 9 Be in- 
duced by the neutron emission into the continuum and 
into bound states of the target nucleus. Both processes 
produce 8 Be which decays into two a-particles, and one 
of them could end up as incomplete fusion. These reac- 
tion mechanisms should be studied theoretically in detail 
in order to develop a better understanding and descrip- 
tion of the breakup, complete and incomplete fusion. In 
fact, a recent experiments 2l| show that the 9 Be breakup 
following neutron transfer dominates the total breakup 
yield. 



VI. CONCLUSIONS 

The complete fusion of 9 Be with spherical target nuclei 
was calculated in the coupled-channels approach. It was 
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FIG. 8: (Color online) Similar to Fig. [7] The calculations 
include, in addition to the decay of the excited states of 9 Be, 
a weak imaginary potential with the strength Wo = -2.5 MeV. 



shown that the B(E2) values for the excitation of the 
ground state rotational band of the 9 Be nucleus can be 
described quite well in the rotor model. This feature 
was exploited in the calculation of matrix elements of 
the interaction between the deformed projectile and a 
spherical target. 

The interaction of the deformed 9 Be projectile with a 
spherical target was calculated using the double-folding 



technique and an effective M3Y NN interaction, which 
was supplemented with an adjustable repulsive term. 
The deformed density of 9 Be was determined so the mea- 
sured quadrupole moment and RMS charge radius were 
reproduced. The densities of the spherical targets were 
calibrated to reproduce the measured charge radii; the ra- 
dius of the neutron density in 208 Pb was calibrated to be 
consistent with the neutron skin thickness predicted by 
Skyrme Hartree-Fock calculations and with the value ex- 
tracted from measurements of antiprotonic 208 Pb atoms. 

The double-folding potential was applied in coupled- 
channels calculations of the fusion of 9 Be with 144 Sm. 
The decay of the excited states of 9 Be was included ex- 
plicitly in terms of complex excitation energies, whereas 
excitations of the target were ignored for simplicity. The 
repulsive part of the effective NN interaction, which es- 
sentially is the only parameter that remains to be deter- 
mined, was adjusted to produce an optimum fit to the 
complete fusion data. The calculated cross sections for 
the decay of the excited states in 9 Be turned out to be 
in very good agreement with the incomplete fusion data. 

Having determined all of the parameters of the theory, 
coupled-channels calculations were performed for the fu- 
sion of 9 Be+ 208 Pb. The complete fusion data were sup- 
pressed at high energies by 20% compared to the pre- 
dicted cross sections, and the incomplete fusion data were 
a factor of three larger than the calculated cross section 
for the decay of the excited states. There are obviously 
other reaction mechanisms, besides the decay of excited 
states, that are responsible for the large incomplete fusion 
cross sections that have been measured for the lead tar- 
get. A likely candidate is the neutron transfer from 9 Be 
to bound states in the target and to continuum states. 
This is also the conclusion of a recent experimental in- 
vestigation by Rafiei et al. [2lj |. 
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